Casimir effect in a wormhole spacetime 



Artem R. Khabibullin a *, Nail R. Khusnutdinov at , Sergey V. Sushkov 6 * 
a Department of Physics, b Department of Mathematics, 
Kazan State Pedagogical University, Mezhlauk 1, Kazan 420021, Russia 
(Dated: February 1, 2008) 

We consider the Casimir effect for quantized massive scalar field with non-conformal coupling £ 
in a spacetime of wormhole whose throat is rounded by a spherical shell. In the framework of zeta- 
regularization approach we calculate a zero point energy of scalar field. We found that depending 
on values of coupling £, a mass of field m, and/or the throat's radius a the Casimir force may be 
both attractive and repulsive, and even equals to zero. 
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O ■ I. INTRODUCTION 

o : 

The central problem of wormhole physics consists of the fact that wormholes are accompanied by unavoidable 
violations of the null energy condition, i.e., the matter threading the wormhole's throat has to possess "exotic" 
properties. The classical matter does satisfy the usual energy conditions, hence wormholes cannot arise as solutions 
of classical relativity and matter. If they exist, they must belong to the realm of semiclassical or perhaps quantum 
■ gravity. In the absence of the complete theory of quantum gravity, the semiclassical approach begins to play the 
most important role for examiningwormholes. Recently the self-consistent wormholes in the semiclassical gravity 
were studied numerically in Refs jlj, [l^, l2(i |23| . It was shown that the semiclassical Einstein equations provide an 
existence of wormholes supported by energy of vacuum fluctuations. However, it should be stressed that a natural 
size of semiclassical vacuum wormholes (say, a radius of wormhole's throat a) should be of Planckian scales or less. 
This fact can be easily argued by simple dimensional considerations [l^. In order to obtain semiclassical wormholes 
having scales larger than Planckian one has to consider either non-vacuum states of quantized fields (say, thermal 
__ 1 states with a temperature T > 0) or a vacuum polarization (the Casimir effect) which may happen due to some 
I/"") I external boundaries (with a typical scale R) existing in a wormhole spacetime. In the both cases there appears an 
<^> additional dimensional macroscopical parameter (say R) which may result in enlargement of wormhole's size. 

In this paper we will study the Casimir effect in a wormhole spacetime. For this aim we will consider a static 
spherically symmetric wormhole joining two different universes (asymptotically flat regions). We will also suppose 
that each universe contains a perfectly conducting spherical shell rounding the throat. These shells will dictate the 
Dirichlet boundary conditions for a physical field and, as the result, produce a vacuum polarization. Note that this 
4^ ■ problem is closely related to the known problem which was investigated by Boyer 8] who studied the Casimir effect 
of a perfectly conducting sphere in Minkowski spacetime (see also |^|). However, there is an essential difference which 
is expressed in different topologies of wormhole and Minkowski spacetimes. A semitransparent sphere as well as 
semitransparent boundary condition were investigated in Refs. 0, El 0> EH Eil • The consideration of the delta- 
like potential which models a semitransparent boundary condition in quantum field theory cause some problems and 
there is ambiguity in renormalization procedure (see the Refs. yUliaHS an d references therein). Thermal corrections 
to the one-loop effective action on singular potential background was considered recently in Ref. p^ |. 

We will adopt a simple geometrical model of wormhole spacetime: the short-throat flat-space wormhole which was 
suggested and exploited in Ref. [20j. The model represents two identical copies of Minkowski spacetime; from each 
copy a spherical region is excised, and then boundaries of those regions are to be identified. The spacetime of the 
model is everywhere flat except a throat, i.e., a two-dimensional singular spherical surface. We will assume that the 
wormhole's throat is rounding by two perfectly conducting spherical shells (in each copy of Minkowski spacetime) 
and calculate the zero-point energy of a massive scalar field on this background. In the end of calculations the radius 
of one sphere will tend to infinity gi ying the Casimir energy for single sphere. For calculations we will use the zeta 
function regularization approach [lfj, which was developed in Refs. [IH SIS G3- 

In framework of this approach, 

the ground state energy of scalar field <j> is given by 



where 



CH S )=E( A (n)+- 2 ) 

(„) 
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is the zeta function of the corresponding Laplace operator. The parameter fi, having the dimension of mass, makes right 
the dimension of regularized energy. The A 2 n j are eigenvalues of the three dimensional Laplace operator C = A — (JZ 

(A - fft)0 (n) = A( n) </>(„), (2) 

where 1Z is the curvature scalar (which is singular in framework of our model, see Eq. ©)• 

The expression Q is divergent in the limit s — > which we are interested in. For renormalization we subtract from 
JU the divergent part of it 

E™ = Vm(E(8)-E**(a)), (3) 

where 

£ div (s) = lim E(s). 

m — >oc 

By virtue of the heat kernel expansion of zeta function is the asymptotic expansion for large mass, the divergent part 
has the following form (in 3 + 1 dimensions) 

IW ^w.-i) (4) 

x jsom 4 r(s - 2) + B 1/2 m 3 T{s - |) + B im 2 T(s - 1) + B 3/2 mT{s - i) + B a r(a)| , 

where -B Q are the heat kernel coefficients of o pera tor C. In the case of singular potential (singular scalar curvature) 
one has to use specific formulae from Refs. [Mil l| for calculation the heat kernel coefficients (see also a recent review 
0). 

Finally, the renormalized ground state energy should obey the normalization condition 

lim E Icn = 0. 

rn — >oo 

For more details of approach see review [|| . 

The organization of the paper is the following. In Sec. [H]we describe a spacetimc of wormholc in the short-throat 
flat-space approximation. In Sec. IHII we analyze the solution of equation of motion for massive scalar field and obtain 
close expression for zero point energy. In Sec. IIVI we discuss obtained results and make some speculations. 

We use units % = c = G = 1. The signature of the spacetime, the sign of the Riemann and Ricci tensors, is the 
same as in the book by Hawking and Ellis [l6[ . 



II. THE GEOMETRY OF THE MODEL 



We will take a metric of static spherically symmetric wormhole in a simple form: 

ds 2 = -dt 2 + dp 2 + r 2 (p)(d6 2 + sin 2 9dtp 2 ), (5) 

where p is a proper radial distance, p G (—00,00). The function r(p) describes the profile of throat. In the paper 
we adopt the model suggested in the Ref. |2J| which was called there as short-throat flat-space approximation. In 
framework of this model the shape function r(p) is 

r (p) = IpI + a > 

with a > 0. r(p) is always positive and has the minimum at p = 0: r(0) = a, where a is a radius of throat. It is 
easy to see that in two regions T> + : p > and 2?_ : p < one can introduce new radial coordinates r± = ±p + a, 
respectively, and rewrite the metric in the usual spherical coordinates: 

ds 2 = -dt 2 + dr 2 ± + r 2 ± {dQ 2 + sin 2 9 dip 2 ), 

This form of the metric explicitly indicates that the regions T> + and D- are flat. However, note that such the change 
of coordinates r± = ±p + a is not global, because it is ill defined at the throat p = 0. Hence, as was expected, the 
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spacetime is curved at the wormhole throat. To illustrate this we calculate the Ricci tensor in the metric iJ^J: 

v a 
n = 2(-l + r' 2 + 2rr") _ J(p) _ 
r 2 a 

The energy-momentum tensor corresponding to this metric has the diagonal form from which we observe that the 
source of this metric possesses the following energy density and pressure: 

-1 + r' 2 +2rr" _ _5(p) 
8nr 2 2na ' 

-1 + r' 2 

r" 8{p) 
Pe = P<p = Z — = 1 — ■ 

III. ZERO POINT ENERGY 

Let us now consider a scalar field <f> in the spacetime with the metric ■ The equation for eigenvalues of operator 
C is 

(A - £72.)^( n ) = A 2 „)0( n ), (7) 

where 7Z is the scalar curvature, £ is an arbitrary coupling with 1Z and A = g Q ^V a V/3, a = 1,2,3. Due to the 
spherical symmetry of spacetime (J5J, a general solution to the equation Q can be found in the following form: 

where Yi n (9, ip) are spherical functions, I = 0, 1, 2, . . . , n = 0, ±1, ±2, . . . , ±Z, and a function obeys the radial 
equation 

u " + 2-u' + ( X 2 — 1 1±}1 -£r)u = 0, (8) 
r \ r 2 I 



where a prime denotes the derivative with respect p, A = y/uj 2 — m 2 and scalar curvature 1Z — —8S(p)/a. For new 
function w = ur this equation reads 

\ r z r / 

and looks like the Schrodinger equation for massive particle with mass M with total energy E = A 2 /2M and potential 
energy 

U = m +-)/2M = 1 —^8{p). (9) 
r am 

Therefore, the £ > 1/4 corresponds to negative potential. 

Unfortunately, in our case it is impossible to find in manifest form the spectrum of operator C given by Eq. JJJ. 
For this reason, we will use an approach developed in Refs. 0, 0, El IE El • This approach does not need an explicit 
form of spectrum. The spectrum of an operator is usually found from some boundary conditions which look like an 
equation W(A) =0 where function \I/ consists of the solutions of Eq. © and depends additionally on other parameters 
of problem. It was shown in Refs. |2tl,'ll ll |6U19| that the zero point energy may be represented in the following form: 

! cos(7ts) f°° j Ufu2 

(n) 



E{s) = _^^ly2 d l ^-^/s-iLin*^), (10) 
2tt 4 J m Ok 
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with the function ^ taken on the imaginary axes. The sum is taken over all numbers of problem and d n is degenerate 
of state [2^. This formula takes into account the possible boundary states, too. If they exist we have to include 
them additively at the beginning in the Eq. (JJJ. But integration over interval |fc| < m (the possible boundary states 
exist in this domain) will cancel this contribution. For this reason the integration in the formula (|10|) is started from 
the energy k — m. Therefore, hereinafter we will consider the solution of the Eq. JSJ for negative energy that is in 
imaginary axes A = ik. The main problem is now reduced to finding the function 'J. Thus, now we need no explicit 
form of spectrum of operator C. 

In the flat regions T>±, where r(p) = ±p + a, r'(p) — ±1, TZ(p) — 0, and in imaginary axes the Eq.© reads 



p ± a 



k 



1(1 + 1) 



(p±a) 



u = 0. 



A general solution of this equation can be written as 

u±[k(a±p)] = A±J^^lAHa±p)] 



2k(a±p) 



K u [k{a±p)}, 



(11) 



(12) 



where /„, K u are the Bessel functions of second kind, v = I + 1/2, and A , B are four arbitrary constants. 

The solutions u^[k(p ± a)] have been obtained in the flat regions T>± separately. To find a solution in the whole 
spacetime we must impose matching conditions for M ± [fc(p± a)] at the throat p = 0. The first condition demands 
that the solution has to be continuous at p = 0. This gives 



u [ka] = u + [ka] . 



(13a) 



To obtain the second condition we integrate Eq.© within the interval (— e, e) and then go to the limit e — > 0. It gives 
the second condition 



du [x] 



dx 



du+\ 



x—ka 



dx 



8i 



u + [ka] 



(13b) 



x—ka 



Therefore, the general solution of Eq. i|ll|) depends on two constants, only. Two other constants may be found from 
Eqs. lfT5a|) and lfl3b)l . 

In addition to two matching conditions i|13aD and l|13b|) we impose two boundary conditions. We round the wormholc 
throat by sphere of radius a + R (p = R) in region V + , and by sphere of radius a + R 1 (p = —R') in region 2?_. 
Therefore the space of wormhole is divided by two spheres to three regions: the space of finite volume between spheres 
and two infinite volume spaces out of spheres. We suppose that the scalar field obeys the Dirichlet boundary condition 
on both of these spheres which means the perfect conductivity of spheres: 



u-[k(R' + a)] = 0, 
u+[k{R + a)} = 0. 



(13c) 
(13d) 



The four conditions (|13|) obtained represent a homogeneous system of linear algebraic equations for four coefficients 
A , B . As is known, such a system has a nontrivial solution if and only if the matrix of coefficients is degenerate. 
Hence we get 



-I v [kc 



-K v [ka] 



l v [ka] 



K v [kc 



KM + ^lAka] Kl[ka] + ^-K v [ka] I' u {ka}- ^I v [ka] K'„[ka] - ^K v [ka] 



I v [k{a + R)] 




K v [k{a + R)] 






I v [k(a + R')} 





K v [k{a + R')] 



(14) 



After some algebra the above formula can be reduced to the following relation for function ^> which we need for 
calculation of the energy iflUjl : 



I v [k(a + R')] (jT 
K v [k{a + R')\ (V 



a- 



K v [ka] 
I v [ka] 



ka 

T 

ka 

T 



K'J\ka 
I' v [ka] 



-K v [k{a + R)] 



(15a) 



-I v [k[a + R)] =0, 



with 



= I v [k{a + R)]K v [ka] - K v [k{a + R)]I v [ka\. 
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In the case R' = R above expression coincides with that obtained in Ref. [2(j. In this case may be represented 
as follows: * m = \f^* 2 , where 

*l = = I v [k{a + R))K v [ka] - K v [k{a + R)]I„[ka], 

*? = U-l^* + ^[I^a + R)}KUka}-KMa + RM[ka}}, 

The solutions of Eq. (|15a|l gives the spectrum of energies between the spheres R and R' . The spectra for regions 
out of these spheres can be found as follows: 

^^K^kia + R)}, (15b) 

* 2 out = K v [k{a + R')]. (15c) 

Indeed, let us consider the energy spectrum of field in space between two spheres with radii R and R > R and Dirichlet 
boundary conditions on them. The solution is a linear combination of two modified Bessel functions 

u r,r = Cil v [kp] + C 2 K„[kp]. 

The Dirichlet boundary conditions give two equations 

dl^kia + R)] + C 2 K u [k{a + R)] = 0, 
Cxl v [k{a + R)] + C 2 K u [k(a + R)} = 0. 

Using these equations we may represent the solution in the following form: 

- - ft ff.|M *»'«' + ft' -*,U. 



™ K,[k{a + R)\ " /„[(.(„ + R) ] 

Let us now assume that R — > oo. In this limit the solution takes the following form: 

uroc = CK v [kp}. 

The Dirichlet boundary condition for this solution on the sphere of radius R gives the equation (|15bl) . As expected 
this condition coincides with expression for space out of sphere of radius a + R in Minkowski spacetime Q . It is 
obviously because the spacetime out of sphere (in general out of throat) is exactly Minkowski spacetime. 
Therefore the regularized total energy i|l(J|) reads 

= a.ooBOra) g y f°° 2 m 2 )V^|- [l n * m + ln^ out + ln*L] . (16) 

I" ^ Jm Ok 

Regrouping terms we can rewrite the above formula in the form having clear physical sense of each term: 

E(s)=AE( S ) + E%(s) + E%(s), (17) 



where 



E%(s) = -^ ^"T ,/ / dk(k 2 ~m 2 )^ 2 - s — ]nl u [k(a + R)]K v [k{a + R)}, (18) 

n , o Jm Ok 



=0 
oo 



E%(s) = -^^^If^v f°° dkie-m^-^lnl^kia + R^K^kia + R')}, 
11 i=q ^ m 

AE(s) = V s ^f, r dk(k 2 ~m 2 )^ 2 - s In* 

f-i Jm Ok 



(19) 
(20) 



and 



I v [k{a + R')]I u [k{a + R)] 



6 



The term E^(s) in the formula (|17f) is nothing but a zero point energy of sphere of radius a + R in Minkowski 
spacetime with Dirichret boundary condition on the sphere 0; note that the term Ej£,(s) has an analogous sense. 

Now we are ready to calculate the Casimir energy for two spherical boundaries by using expression Qlfijl and Eq. 
©. Then let us consider the Boyer's problem. We consider "gedanken experiment": we take a single conducting 
sphere and measure the Casimir force in this situation. For this reason we have to take a limit R' — > oo. In this case 
the energy (|19J) tends to zero, and so the term AE(s) in Eq. (|2C)|> represents the difference between Casimir energies 
of a sphere rounding the wormhole and a sphere of the same radius in Minkowski spacetime without wormhole. In 
the limit R' — > oo we find 

* - - '^m§) (H) ^ + t*m) - im§ <»> 

If one turns R — > oo then the energy E^ tends to zero and so 



KuM (U~l) K v [ka] + !jK' v [ka]\ . (22) 



This expression coincides exactly with that obtained in Ref. |20j| and describes the zero point energy for whole 
wormhole spacetime without any additional spherical shells. 

A comment is in order. As already noted the positive £ corresponds to attractive potential and therefore the 
boundary states may appear. The appearance of boundary states with delta-like potential has been observed in Ref. 
[2l| . Thus, we have to take into account the boundary states at the beginning. Nevertheless, the final formula i|16|) 
contains these boundary states, as it was noted in Ref. But it is necessary to note, that in this paper we will 
consider £ < 1/4. Indeed, let us consider for example I = 0. In this case 



«I» = - e ~ ka cosh(fc(a + R)) { cosh(kR) + 



2^ 



sinh(fcii) 



For £ > 1/4 this expression may be equal to zero for some value of k > m, R and a and integral i|16|) will be divergent. 
As noted in Ref. |2l| in this case we can not use the present theory. The same boundary for £ was noted in Ref. 
[20| . This statement is easy to see from expression for potential energy given by Eq. JJjJ . For £ > 1/4 the energy is 
negative and the boundary states may appear. 

The general strategy of the subsequent calculations is following (for more details see Refs. 0, IS El IS H])- To single 
out in manifest form the divergent part of regularized energy we subtract from and add to integrand in Eq. (|16[) its 
uniform expansion over 1/v. It is obviously that it is enough to subtract expansion up to 1/v 2 , the next term will 
give the converge series. We may set s = in the part from which we had subtracted the uniform expansion because 
it is now finite (see Eq. The divergent singled out part will contain the standard divergent terms given by Eq. 

PJ) and some finite terms which we calculate in manifest form (all terms except A in (J23J) ) ■ 

The uniform asymptotic expansions both 121|) and 122fl are the same for R ^ 0. Indeed, in this case the ratios 



I„[ka\ K v [k{a + R)} _ 

2u ln(l+ — } 

K v \ka\ I u [k{a + R)] 

1 K v [k{a + R)} _ 2l , ln(1+ H) 
Kl\ka] I v [k{a + R)] ~ 

are exponentially small and we may neglect them. The well-known uniform expansions of Bessel functions [l| were 
used in these expressions. For this reason we may disregard this fraction in Eq. (|21|l and arrive to Eq. I|22l) . This 
is a key observation for next calculations. Due to this observation the divergent part which we have to subtract for 
renormalization from l|2Q(l has been already calculated in Ref. 20] . By using the results of this paper we may write 
out the expression for rcnormalized zero point energy: 

AE = -^L_(Mn^ + 0), (23) 

3 

il = A+ MP), (24) 
fe=-i 

b = \b^ -h(3 2 +b 2 , (25) 
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FIG. 1: The plots of renormalized zero-point energy E re n/m as a function of x = R/a for (5 — 0.04,0.16,0.5 and for various 
values of £ and fixed mass m. We observe that increasing £ leads to appearance maximum and/or minimum. For subsequent 
increasing £ the curve will turn over and extremum disappears. If the radius of spherical shell exceeds ten radius of throat the 
zero-point energy takes on a value which equals to zero-point energy in whole wormhole spacetime. 



1=0 h 



where 

dy\ y 2 - fin* + 2vr,{y) + -N ± - \n 2 + ^N 3 ) , (26) 

p/u v v oy \ v v 2 v 6 ) 

* = l^ N] -^ N] /,N(i + ,)]J sj + T^^J - 8 + ■ (27) 

&fe are the heat kernel coefficients, (3 = ma is a dimensionless parameter of mass, and x — R/a is a dimensionless 
parameter of sphere's radius. The explicit form of heat kernel coefficients b^, and also expressions for lu^, N^, r\ 
are given in the Ref. j2£j. Note that they do not depend on the radius of sphere R. The only dependence on R is 
contained in the coefficient A which has to be calculated numerically. The expression for contribution of the sphere 
in Minkowski spacetime l|18[) may be found in Ref. We only have to make a change R — > a + R. 



IV. DISCUSSION AND CONCLUSION 



In this section we will discuss results of numerical calculations of zero-point energy given by formula 1)23(1 . The 
renormalized zero-point energy is represented in figures QEl as a function of x — R/a (the position of sphere rounding 
the wormhole) for various values of (3 = ma and £. (Note that the value x — R/a characterizes the position of sphere 
rounding the wormhole; x — corresponds to sphere's radius equals to throat's radius.) In Fig. ^we only show the 
full energy E. Note that the AE differs just slightly from the full energy E. For the same reason we reproduce in 
Fig. |3 the AE, only. 

Characterizing the result of calculations we should first of all stress that the value of zero point energy E ren in the 
limit R — > oo tends to some constant value obtained in Ref. 20] for the case of wormhole spacetime without any 
spherical shells. In the limit R — > (i.e., when the sphere radius a + R tends to the throat's radius a) the zero-point 
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FIG. 2: The plot of AE ren /m as a function of x — R/a for £ = ~ and for various values of /3 and fixed mass m. We observe 
the dynamics of deformation of energy due to changing the parameter (3 = ma for fixed £. 

energy E ren is infinitely decreasing for all /3 and £. This means that the Casimir force acting on the spherical shell 
and corresponding to the Casimir zero point energy E ren is "attractive", i.e., it is directed inward to the wormhole's 
throat, for sufficiently small values of R. In the interval < R/a < oo there arc three qualitatively different cases of 
behavior of E ren depending on values of (3 and £. Namely, (i) the zero point energy E ren is monotonically increasing 
in the whole interval < R/a < oo. There are neither maxima no minima in this case. Hence the Casimir force is 
attractive for all positions of the spherical shell, (ii) E ren is first increasing and then decreasing. A graph of the zero 
point energy has the form of barrier with some maximal value of E ren at R\/ a. The Casimir force is attractive for the 
sphere's radius R < R\ and repulsive for R > R\. The value R = R± corresponds to the point of unstable equilibrium, 
(iii) The zero point energy E ren is increasing for R/a < Ri/a, decreasing for R%/a < R/a < R2/a and then finally 
increasing for R/a > Rz/a, so that a graph of E ren has a maximum and minimum. In this case the Casimir force is 
directed outward provided the sphere's radius R\ < R < i?2, and inward provided R < R\ or R > i?2- Now the value 
R = i?2 corresponds to the point of stable equilibrium, since the zero point energy E ren has here a local minimum. 

It is worth noting that the Casimir force is attractive in the whole interval < R/a < oo for sufficiently small 
values of £ and/or large values of p. Otherwise, it can be both attractive and repulsive depending on a radius of 
sphere rounding the wormhole's throat. The similar situation appears for delta-like potential on the spherical or on 
the cylindrical boundaries [25L |2(| . The repulsive Casimir force was also observed in Ref. ^3 f° r scalar field living in 
the Einstein Static Universe. 

The considered model let us speculate in spirit of Casimir idea who suggested a model of electron as a charged 
spherical shell |!j ■ Casimir assumed that such a configuration should be stable due to equilibrium between the repulsive 
Coulomb force and the attractive Casimir force. However, as is known, this idea does not work in Minkowski spacetime 
since the Casimir force for sphere turns out to be repulsive @- Now one can revive the Casimir's idea by considering a 
spherical shell rounding the wormhole. In this paper we have shown that the Casimir force now can be both attractive 
and repulsive. Moreover, there exists stable configurations for which the Casimir force equals to zero; the radius of 
spherical shell in this case depends on the throat's radius a as well as the field's mass m and coupling constant £. 
Thus, one may try to realize the Casimir's idea taking a sphere rounding a wormhole. Of course, our consideration was 
based on the very simple model of wormhole spacetime. However, we believe that main features of above consideration 
remain the same for more realistic models. 
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